Radiation pressure driven vibrational modes in ultra-high- Q silica microspheres 
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Quantitative measurements of the vibrational eigenmodes in ultra-high-Q silica microspheres are 
reported. The modes are efficiently excited via radiation-pressure induced dynamical back-action of 
light confined in the optical whispering-gallery modes of the microspheres (i.e. via the parametric 
oscillation instability) . Two families of modes are studied and their frequency dependence on sphere 
size investigated. The measured frequencies are in good agreement both with Lamb's theory and 
numerical finite element simulation and are found to be proportional to the sphere's inverse diameter. 



Silica microcavities [J such as microspheres [2] or 
microtoroids[3] possess ultra-high-Q optical whispering 
gallery modes (WGMs), while simultaneously exhibit- 
ing mechanical modes which lie typically in the radio 
frequency range. Owing to the resonant buildup of 
light within these cavities the effect of radiation pres- 
sure is enhanced, leading to mutual coupling between 
the mechanical and optical modes, as first predicted by 
Braginsky in the context of the Laser Interferometer 
Gravitational Wave Observatory, LIG0[4]. When en- 
tering the regime where the photon lifetime is compa- 
rable to the mechanical oscillation period and the cav- 
ity is pumped with a laser whose frequency slightly ex- 
ceeds the WGM resonance (i.e. blue-detuned excitation), 
this mutual coupling gives rise to a parametric oscillation 
instability!!] which is characterized by regenerative me- 
chanical oscillation of the mechanical eigenmodes. This 
phenomenon has been first reported in toroid microcav- 
ities [51 |5J [7j. On the other hand, red-detuned light 
can induce cooling of the mechanical mode, as recently 
reported[8l [9l [10]. In this letter, parametric oscillation 
instability in ultra-high-Q silica microspheres is observed 
and the mechanical resonant frequencies and mode pat- 
terns studied. In contrast to earlier studies of acoustic 
modes of nanospheres[TTJ [12] using Raman or Brillouin 
scattering from ensembles, the present method allows 
measurement of the mechanical modes of single micro- 
spheres in a larger diameter regime (35-110 /im in our 
case). Furthermore, the mechanical Q-factors are deter- 
mined. 

Ultra high-Q {Q > 10^) silica microspheres are fabri- 
cated by melting the tip of a single mode optical fiber 
with a CO2 laser (A = 10.6 /zm). Due to the strong ab- 
sorption of silica around 10.6 ^m, surface tension induces 
highly symmetric silica spheres with a near- atomic ally 
smooth surface [1 3 . The sphere is held by a thin fiber 
stem (cf. Fig. 1 inset). The WGM are excited with high 
ideality ^14, by evanescent coupling via a tapered optical 
fiber [IS] using a 1550— nm tunable external cavity diode 
laser as pump source. Owing to the high finesse (ap- 
proximately 10^), the large optical energy stored in the 
microcavity exerts a force on the cavity sidewalls due to 
radiation pressure. This force can give rise to regenera- 
tively driven mechanical oscillations if photon lifetime is 
similar to the inverse acoustic resonance frequency [3 [6]. 
In essence, the radial force exerted by radiation pres- 
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FIG. 1: Spectrum of the transmitted optical laser of the fiber 
taper when coupled to a silica microsphere with a diameter of 
49 /im. Two families of spheroidal mechanical modes could 
be driven regeneratively. Specifically, the modes are indenti- 
fied as i^(i,2) (a quadrupole mode) and 1^(10) (a radial breath- 
ing mode). The launched power in this experiment was 600 
micro-Watts and was sufficient to exceed the threshold for 
parametric oscillation instability. The inset shows the SEM 
image of the microsphere. 



sure takes the cavity out of resonance by deformation 
of the cavity wall which causes subsequently a reduction 
in the radiation pressure force. The whole process re- 
sumes upon the restoration of the original shape of the 
cavity, leading to a periodic motion of the cavity. In 
this way, modes with radial deformation can be excited 
which modify the path length of the optical wave and 
therefore affect the magnitude of the radiation force. In 
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FIG. 2: Experimentally measured frequencies of the 
spheroidal modes, with stars denoting the !^(i,2) mode and 
dots denoting the i^(i,o) mode. Numerically calculated eigen- 
frequencies of these modes are shown in blue solid line (i^(i,2)) 
and red dashed line (2^(1,0) )• Spheroidal mode with the 

frequencies lie between /^(i,2) and !^(i,o) j which is not ex- 
perimentally observed, is also presented in black dotted line. 
The inset shows the relationship between the frequencies of 
the eigenmodes and the inverse diameter of the silica micro- 
spheres. 



the present experiments the threshold for the parametric 
instability was in the range of typically 100 micro- Watts. 
The driven mechanical oscillations causes the appearance 
of motional sidebands (and their harmonics) which can 
be readily detected in the spectrum of the transmitted 
laser light using an electronic spectrum analyzer. Fig. 1. 
shows a typical spectrum of the transmission, showing re- 
generative oscillations of two different mechanical eigen- 
modes of the same sphere under different taper loading 
conditions. By changing the taper loading (and hence 
the optical Q) the value of the oscillation threshold of 
the two modes crosses, thereby causes a switching from 
the low- to a high frequency mode as coupling strength 
increases and the optical linewidth decreases, in agree- 
ment with the theoretical predictions [IB]. To identify 
the mode families, the vibrational frequencies of the two 
lowest lying mechanical frequencies were recorded as a 
function of size. While observable in principle, light- 
induced modifications of the mechanical modes' dynam- 
ical properties lOj, in particular shifts in the mechanical 
resonance frequency, are considered to be small (relative 
frequency shift typically < 0.1%) and are neglected in 
this letter. The result of this study is shown in Fig. 2. 
As evident, the spheres mechanical frequencies are in- 
versely proportional to the sphere diameter. 

Next, the observed modes where identified by numer- 
ical studies. Since the fabricated silica microspheres ex- 
hibited no observable eccentricity under SEM imaging 
(see the inset of Fig. 1), and the diameter ratio of the 
stem holding the sphere and the sphere itself is on the 



order of 0.1, the sphere can be considered as almost free; 
hence it is a judicious choice to adopt the stress-free 
boundary condition. Studies on the nature of the funda- 
mental modes of vibration for small elastic spheres with 
free-surface boundary condition are well known. The first 
well-established theory was formulated by Lamb, with 
two types of modes predicted, the spheroidal and tor- 
sional modes [17j. The equation describing the wave prop- 
agation in a homogeneous elastic body with free surface 
can be written as [181, 



pu= {X + 2^)V(V • u) - X (V X u) 



(1) 



where u is the displacement vector, p is the mass 
density, and A and fi are Lame constants. Here A = 
(i+a)ii^2a) ' M = 2{i+a) ' "^i*^ ^ denoting the 
Young's modulus and a the Poisson ratio of the mate- 
rial. Eq. (1) can be solved by introducing a scalar po- 
tential (j>o and two vector potentials $1 = {rcfii, 0,0) and 
$2 = {r(f>2,0, 0) with u = V(/)o + Vx$i-|-VxVx$2- 
Then the general solutions of the equations resulting from 
(1) are written as 
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where ji is the spherical Bessel function and F/" is the 
spherical harmonic function and Vq is the longitudinal 
sound velocity and Vi = V2 are the transverse sound 
velocities. An angular momentum mode number £ ( 
^ = 0,1, 2..), an azimuthal mode number m (—£ < m < £) 
and a radial mode number n (n — 1,2,..) are used to 
characterize the acoustic modes, where n = 1 corre- 
sponds to the surface mode and n > 2 to inner modes 
and i'n,i.m denotes the frequency of the vibration charac- 
terized by the mode numbers (n,£,m). It is noteworthy 
that a spheroidal mode with angular momentum £, is 
{2£ + l)-fold degenerate, hence in Fig. 1 we use {n,£) 
instead of (n, £, m) to assign the eigenfrequencies. 

Two classes of modes are derived when applying the 
free boundary condition|18j. One of them is the torsional 
vibration which induces only shear stress without volume 
change, and no radial displacement takes place in these 
modes. Thus these modes cannot be excited using ra- 
diation pressure which relies on a change in the optical 
path length. In contrast, the class of mode in which vol- 
ume change is present is referred to as spheroidal modes. 
According to Lambs theory the £ = spheroidal mode 
eigenvalue equation is written as 



tan(hR) 
hR 



1 - Uk'^/h^)h'^R^ 







(3) 



where k — 27ri//Vi, h — 2TTiy/Vo, R is the radius of 
the sphere, Vq — ^(A -f 2^)/p, and Vi — V2 — \/JiJp. 
Other eigenvalue equations for torsional modes and ^ > 
spheroidal modes are contained in Ref.jT?]. Next, the res- 
onant frequencies of the lowest frequencies (i.e. n = I, 
£ = 0, 1, 2) were numerically calculated as a function of 
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FIG. 3: Finite element modeling of three spheroidal modes 
i^(i,2) (left) and z/{i,o) (right) of a silica microsphere with its 
von Mises stress (color coded) and deformed shape (greatly 
exaggerated for clarity). 



sphere size (compare Fig. 2, solid lines). The eigenval- 
ues versus the inverse of the diameters are shown in the 
inset. As seen from Fig. 2 the measured data fit very 
well to the theoretical prediction based on the V(ifi) (ra- 
dial breathing) and J^(i,2) (quadrupole) mode, and reveals 
that the eigen- frequencies of the microspheres have linear 
dependence on the inverse microsphere diameter. 

We note that the sphere's eccentricity can lift the 
{2£ + l)-fold degeneracy in the mode number m of the 
^{1,2) mode[19j, which in the experiments was not ob- 
served owing to the high degree of symmetry of the mi- 
crospheres. 

It is worth pointing out that £ — 1 spheroidal mode 
has the eigenfrequency lying between i — and £ — 2 
spheroidal modes (Fig. 2); however, it is not observed 
in the experiment since £ = I spheroidal mode result in 
a vanishing path length change around the optical whis- 
pering gallery mode trajectory such that the mutual cou- 
pling of optical and mechanical mode vanishes. 

To gain a more complete understanding of the mechan- 
ical modes, the numerical studies were complemented 
with finite element simulations of the stress and strain 
fields. Using axial symmetric finite element modeling, 
the mode families ^^(1,0) ^^'^ ^^(1,2) were calculated and 
excellent agreement found with the numerical solution 
of the preceding section. In the simulation, a Young's 
modulus of 73.1 x 10^ Pa, a Poisson ratio of 0.17 and a 
density of the microsphere of 2.203 x 10'^ kg/m? is used 
(Corresponding to A = 16.09x10^ Pa and /i = 31.24x10^ 
Pa) . The longitudinal and transverse sound velocities are 
5972 m/s and 3755 m/s, respectively. Fig. 3 depicts the 
deformed shape and von Mises stress, as obtained from 
finite element simulations of two spheroidal modes with 
the lowest frequencies. The stresses are greatly exagger- 
ated for clarity. 

To determine the mechanical dissipation of the vibra- 
tional modes, the mechanical quality factor of the sphere 
modes where measured in the sub-threshold regime (in a 
purged nitrogen environment). For this measurement the 
laser power was adjusted to a level far below the threshold 
of the mechanical oscillation. In addition, the mechan- 



ical Q factor was measured for both blue-detuning (i.e. 
where the radiation pressure decreases the mechanical 
dissipation, causing mechanical amplification) and red- 
detuning (where the radiation pressure causes the me- 
chanical damping to be increased, leading to cooling). 
Both values coincided closely, and yielded Q- values of up 
to 7000 for the radial breathing mode. Note that the Q 
was difficult to measure for the £ = 2 mode owing to the 
5-fold degeneracy of the mode which precludes fitting of 
the spectrum with a single oscillator function. 

In conclusion, spheroidal acoustic modes in silica mi- 
crospheres driven by radiation pressure induced para- 
metric oscillation instability are reported. The observed 
vibrational modes are indentified as spheroidal modes 
whose frequencies agree well with Lamb's theory and nu- 
merical simulation and exhibit a linear dependence on 
the inverse of the sphere diameters. 
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